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The circuit complexity of quantum qubit system evolution as a primitive problem in quantum computation 
has been discussed widely. We investigate this problem in terms of qudit system. Using the Riemannian 
geometry the optimal quantum circuits are equivalent to the geodetic evolutions in specially curved 
parametrization of 5[7(<f°). And the quantum circuit complexity is explicitly dependent of controllable 
approximation error bound. 



Quantum algorithms based on the circuit model have been explored to solve difficult problems in terms of 
classical computers', e.g., large integer factoring^ ' and quantum searching algorithm"*. These quantum 
models are composed of a sequence of quantum gates which may be generated from unitary evolution 
operator U{t) with special Hamiltonian. The interacted gates or nonlocal quantum operations are views as the key 
of circuit complexity. One algorithm is computable or efficient if it can be realized with polynomial interacted 
gates in terms of system size^'"". One main goal of quantum computation is to investigate efficient quantum circuits 
to synthesize quantum operation^"^. 

The quantum circuit complexity for qubit systems has been investigated with different approaches such as the 
geometric approach'" It provides an equivalent statement as the shortest path problem in a curved geometry. 
The quantum circuit complexity may be reduced to the metric distance between the identity operation and 
desired unitary using designed Riemannian metric on the unitary space. Thus the metric distance is a good 
measure of the complexity of synthesizing quantum nonlocal evolution. 

In comparison to the qubit system, d-dimensional quantum states (qudits for short) could be more efficient in 
quantum applications. With larger state space the qudit algorithms may improve channel capacity'^"" and 
quantum gates implementation" '^, increase security^""^'' and explore quantum features^^"^^. The qudit systems 
have also been experimentally realized^"". The high-dimensional quantum system maybe provide different 
quantum correlations and efficient information processing, d degenerate (in the rotating- wave approximation 
sense) ground states may form the qudit, coherently coupled via a common excited state by pulsed external fields 
of the same time dependence and the same detuning^^. The multiple photon states may also be reformed as qudit 
with the help of symmetric primitive states under the permutation invariance''^ '^. 

In this paper we focus on the circuit complexity of qudit quantum computations. The results for qubit- 
systems'-' and qutrit-system'* are extended to general qudits. Algebraic methods such as eigenvalue decomposi- 
tion''* ''^, cosine-sine decomposition"" or Householder Reflection decomposition''^ and qudit computation based 
on d-level cluster^" have been used to find the optimal qudit circuits. In comparison to these evaluations, the 
quantum circuit complexity of nonlocal qudit evolutions are revalued in terms of the Riemannian geometry on 
SU{d"). The qudit system evolution is completely determined by the geodesic equation and the initial value of the 
Hamiltonian. The minimal geodesies through a desired unitary, will be investigated by deforming the known 
geodesies of a special metric to the geodesies of the metric of qudit circuits. Thus the qubit or qutrit systems"''' are 
improved to qudit ones, in the sense that universal qudit gates are possible to efficiently synthesize a unitary 
operation with controllable error, and present explicit relationship between the quantum circuit complexity and 
total error bound, which have not ensured in previous results'''''''. Three-qudit system as an explicit example is 
presented in supplementary information. 
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Results 

A quantum gate on «-qudit states is a unitary matrix U e SU{d") 
determined by time-dependent Hamiltonian H{t) (with reduced 
Planck constant) according to the following Schrodinger equation, 



i-U{t)=H{t)U{t) 



(1) 



with 17(0) = 7, [/(T) = 17 and ! = V - 1 . The evolution time t ^ [0,T]. 
H can be decomposed with the generalized GeU-Mann matrices. The 
Lie algebra su{d") of SU(d") is different from the Pauli matrices of 
qubits. Thus we first present necessary results of SM(d"). 

Let fiji- denotes the d-hj-d matrix with a 1 in the j, k position and 
zeros elsewhere, a basis can be described as follows: 



--ejk + etj,! <j<k<d, 



Ujk = i[ejk -ekj),l<k<j<d, 



M^l = diag(l, ■ ■ ■ ,1, -;,0rf_2j),i= 1, ■ ■ ■ ,d- 1 



(2) 



(3) 



(4) 



Here, diag represents the diagonal matrix, 0^; _ 2/ denotes the zeros of 
length d ~ 2j. Thus traceless and Hermitian matrices 

which are generalizations of Pauli matrices for qudit systems. 
These matrices combined with identity matrix : = u'^^^ span the 
vector space of d X d Hermitian matrix 

/ Xu Z12 Zl3 ■ ■ ■ Zld \ 
Zu X22 Z23 ■ ■ ■ Z2d 

Zn Z23 ^33 ■ ■ ■ Zid 



\ Zld Z2d Zid 



(5) 



Xdd / 



with real Xjj and complex Zjk, and are also named with generalized 
GeU-Mann matrices. They construct one group representation of 
su{d) and other representations may be obtained by arbitrary unitary 
transformations. From the unitary matrix representation a recursive 
parametrization of the unitary matrix Ud is followed as 



Ud-- 

with special coefficients Hji^. 



Let X( = M -J. with £ =jd + k and 



Xj = I®'-i(x)xf(g)7®''^ 



(6) 



(7) 



be an operator acting on the s-th qudit with Xf and the rest qudits 
with identity operation I, £=l, - ■ ■ . The basis of su{d") is consti- 
tuted by {yf'}, t=l,---,n, Pt = {ii, ■ ■ ■ ,it} with all possible 
1 < ii < ■ ■ ■ <ik<n, where 



k=i 



(8) 



denotes all operators with generalized GeU-Mann matrices 
Xj^,---,Xji^ acting on t qudits at sites P = {ii, ■ ■ ■ ,ik\, respectively, 
and rest with identity. The element in { Y^' }• may be named as (-body. 

To evaluate the error bound"''"* the norm of an operator A is 
defined by 

||A||= maxllAxIl (9) 
lkll = i 

which is equivalent to the operator norm given by (A, B) = tiA' B. 
The norm of generalized GeU-Mann matrices satisfies 



|uj)t|| = l,j7^A:, and ||mj:(:||=A: with fc=l,- 
1 



,d—l. If we replace 



these not normalized lij-j- with - ui^i^, the generalized GeU-Mann mat- 

K 

rices are then normal with respect to the operator norm'"*, stUl 
denoted by { Y^' }• , is normalized. 

Our main result is the following theorem 

Main theorem. For any small constant e, each unitary Ua g SU(d") 
may be synthesized using 0(tr^j one- and two-qudit gates, with error 
\\U-UA\\<e.. 

To complete the proof necessary lemmas should be proved. Notice 
that the generalized GeU-Mann matrices have the following com- 
munication relations [xj, x^] = 0 or c^-pCj for suitable constants Cji^ 
and generalized GeU-Mann matrix x^, speciaUy, each generalized 
GeU-Mann matrix x^ may be generated by the Lie product of two 
other generalized GeU-Mann matrices Xj and Xj^. By using the tensor 
product (xi ® X2){x3 ® X4) = (X1X3) ® (X2X4), it easily follows that 
Xi, ®x,®xe^ = [xt^ ®{l/y^)xj®IdJd® {'i-/y^)xk®Xi^] . Thus 
all 3-body interaction may be generated by 1-body and 2-body. By 
induction it is not difficult to prove the following lemma. 

Lemma 1. The 1-body and 2-body are universal for SU{d"). 

This lemma may be reformed as: aU s-body items (s s 3) in the 
basis of su{d") can be generated by the Lie bracket products of 1 -body 
and 2-body items. With the lemma only problem is to address the 
approximate accuracy under the limited 2-body items. 

The time-dependent HamUtonian H{t) is represented as 



H=Y.h,aj+Y,hsCT[ 

jeI, SEI2 



(10) 



where Ii denotes all possible one and two-body interactions, I2 
denotes all other more-body interactions, hj are real coefficients, 
i.e., cij denotes aU possible Ff ' , yf^ defined in equation (8) and a[ 
denotes all possible Yf'{t>3i) defined in equation (8). The cost of 
unitary operation U synthesization with Hamiltonian is defined as, 
simUar to the qubit case. 



F{H)^(j2h]+pJ2hl 

\jeli still 



1/2 



(11) 



where p > 0 is the penalty paid for three-and more-body items"'". 
This cost gives rise to a natural distance in SU{d"). A curve Cu 
between the identity operation / and L7 is a smooth function, 

r 17(0 : [0,T]^SU{d") 



1 U{0)=I,U{T) = U 



(12) 



The curve length is given by D(Cu) = j F{H{t))dt. D{Cu) is invari- 

J 0 

ant under different parameterizations of Cu, letting F{H{t)) = 1 by 
rescaling H{t) the evolution time T — D{Cu)- The distance D{I, U) 
between / and U is defined by 

(13) 



D(I,Lf)= min d{Cu) 

vCii 



The function F{H{ • )) has defined a norm associated to a right invari- 
ant Riemannian metric'^' whose metric tensor 0 is presented as 

' 0, if (Ti 7^ (72 

0tri,ff2 ~ > 1' "^1 ~'^2 and o", are 1 or 2 — body (14) 
p, if ffi = o"2 and o", are no less than 3 — body 

With this metric the distance D{I, U) equals to the shortest path of the 



geodesic equation (^—H{t)J{t)J=i{H{t),[H{t)J{t)]). Here <•, •) 
is the inner product on the tangent space su{d"), and /is an arbitrary 
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operator in su{d"). The Lemma 1 with the basis { Y'^' } of su(d") shows 
all q-body items (q > 3) may be generated by Lie bracket products of 
1-body and 2-body items. Thus the metric in equation (14) is reas- 
onable to find the minimal length solution to the geodesic equation'^', 
because the multiple-body items will be ignorable for large p. The 
one- and two-body items mainly contribute to the solution. 

Based on the statement above the HamUtonian H{t) may be pro- 
jected onto one- and two-qudit items Hjit). The projection error is 
controllable in terms of large penalty p. And then the evolution 
according to H2{t) is divided into small time intervals and approxi- 
mated by a constant mean HamUtonian over each interval. 
Moreover, the evolution according to the constant mean 
HamUtonian over each interval may be composed of a sequence of 
one- and two-qudit quantum gates'"* '*. The main problem is to evalu- 
ate the total errors introduced by these approximations. 

Lemma l.IfU ^ SU(d") is generated by H(t) satisfying \\H(t)\\rSc in 
time interval [0, T], then 



(15) 



where H = 



T 



H(t)dt is the mean HamUtonian. 



Proof. The unitary operator U{t) has defined one Dyson operator^"* 
with the evolution = U{t)\>j/(Si)). Its Tomonaga-Schrodinger 

equation is defined 

■'^■u{t)\m)=mm\m) (i6) 



dt 



U. This 



with U{t)=I-i[ H{ti)U{ti)dti and U{T) = 
J 0 

representation leads to the Neumann series of U(T) as 

U{T)=f2(- '■)" f t tl)H(t2) • • • H{t„)dt,dt2 ■ ■ • A:( 17) 

m=a J 0 J 0 Jo 



with ti>t2> 



Choosing ti< ■ 



>t„,. The second term is 
T 



-II H{ti)dti = -iHT. 

0 



for 



an mteger 



T fh 



oJo Jo 

m— 1\ 



dtidt2 ■■■dt„ = T"/{m\L'") because of mj! 



then 



r"'/(m!r 

\\e-"''-U 



Thus we have 



■dt„ 



(fT" , Lc"'T"'\ (18) 



= e''^ -l-cT + Le''^/'--L-cT 
for L ^ CO, from the standard norm inequality | \XY\ \ 

mm - c. 



||X||||Y|| and the condition 



Lemma 3. Suppose U2 s SU(d") is a unitary matrix generated by 
H2(t) with evolution time T. Then 



\\U-U2\\<d"p-'/'DiCu) 



(19) 



Proof. Assume two time-dependent HamUtonians H{t) and H{t) 
(with reduced Planck constant) generate U{t) and V{t) respectively 



i- U{t)=H{t)U{t),i-V{t) =H{t)V{t) (20) 

with L7(0) = V(0) = I, U{T) = U and V{T) = V. Using the triangle 
inequality and unitary invariance of operator norm, we obtain 

\\u-v\\ = \\v*u-i\\ 



= f\\iV{tr{H{t)-H{t))U{t)\\dt 
J 0 



f \\V{tr{Hit)-H{t))U{t)\\dt 
J 0 



(21) 



\H{t)-H{t)\\dt 



Here, l7-y= f i{H{t)V{t) -H{t)U{t))dt is used. The Euclidean 
■I 0 

norm of H{t) = hjOj is given by ||//(t) = _ j h]- From the 
Cauchy-Schwartz inequality it follows 



iimoii = 



1/2 

<EH<d"{J2hj] =d"\\H\\ 



(22) 



Therefore from equations (21) and (22) the distance 
D{Cu)= rp{H{t))dt 



> ( F[H{t)-H2{t))dt 
J 0 

.J 0 

>'^j\H{t)-H2{t)\\dt 



(23) 



Lemma 4. The distance D(I, U) always has a supremum independent 
ofp in terms of large p. 

To prove this lemma, the Chow's theorem'"'' is required. Let M be a 
connected manifold and Q a connection on a principal G-bundle. The 
Chow's Theorem states that the tangent space = HyM @ VJSA at 
any point q e M with the horizontal space Hjvl and vertical space 
VJVl {VxM=g, g denotes Lie algebra of G). Denote {Z^} as a local 
frame of H^M, then any two points on M can be connected by a 
horizontal curve if iterated Lie brackets of {Zj"} evaluated at q span 
the tangent space M^^. 

Proof of lemma 4. The connection and completeness of SU(d") gives 
one choice to find the tangent space Mj as its Lie algebra su{d") at the 
identity matrix /. Thus I and U e SU{d") is connected by a unique 
geodesic link associated the Riemannian metric in equation (14). The 
distance D{I, IT) is monotonically increased in the penalty p. 
Moreover, from the Lemma 1 su{d") may be generated with 
iterated Lie brackets of 1-body and 2-body items in the basis 
{Y^'}. With the Chow's Theorem the curve connecting / and Uin 
horizontal subspace is unique in terms of the metric in equation (14) 
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because of the invariance of su{d"). Otherwise, special geodesic with 
initial tangent vector in su(d") may be existed. Hence D{I, U) has an 
optimal upper bound do independent of p. 



Lemma S". If A and B are two unitary matrices, then 

||A*-B*||<fc||A-B|| 



(24) 



for any integer k. 

Based the unitary invariance this lemma is easily followed from 



<||A-5|1 
<\\A-B\\ 
= k\\A-B 



|a'''"M| + ||a*"^b|| + 



+ ||B*-^||) 



(25) 



Lemma 6. Suppose H(t) only contains one- and two-body 
Hamiltonian with \hj\ si at time interval [0, T]. Then there is a 
unitary Ua satisfying 

\\e-'"'^ -UA\\<cn'T'-N-^ (26) 

and being synthesized using at most Nd*n^ one- and two-qudit gates, 
where c is a constant and N is large integer. 

Proof. Divide the interval [0, T] into JV intervals of size T/N. In every 
interval, define a unitary matrix 



u- 



T/N~- 



27 



where! = d^{d^ — l)(n^ — w)/2 = 0{n^) (tf — 1 generators ofsM(d")) 
denotes the total number of possible one- and two-body interactions 
in Hit). From the Trotter formula*" (e"''* + ^" = e""e"" + 0{t^) for two 
Hermitian matrices A and B) there exists a constant c such that 



-'{HU'''''' 



T/N\ \ — 

<cn^T^N-^ 



(28 



Using the Lemma 5 

1 1 g-'-Hi- _ ^ 1 1 < n| I _ Ut,^ 1 1 < r^N- 1 (29) 

which means that one can approximate e'"^ with at most NL £ Nd*n^ 
quantum gates. The error is controlled by the division number N. 
With these lemmas we can prove main Theorem. 

Discussion 

The present Theorem have explained the circuits complexity of in 
quantum computation with n-qudit systems in terms of the 
Riemannian geometry. Similar to the qubit case'^ " and qutrit case'"* 
the optimal quantum circuit is reduced to the shortest path problem 
based on special curved geometry of SU{d"). The qudit systems pre- 
sent different algebraic derivations from qubits and qutrits. 
Especially to realize the unitary invariance of the norm we take use 
of operator norm in equation (9) similar to the qutrit case while the 
norm ||M||i = max^^|^) = ilK'^j-Mji^)!} used before is not unitary 
invariant'^'". For instance, consider 



(30) 



O denotes zero matrix. One has ||M||i = 1/2 and ||MLf||i = 
imXx2+y:i\axy + b*y^\ ^ 1/2. Generally, one has j|M||i s jjM|| from 





/O 1 






/a b 


M = 


0 0 


J 




b* —a* 










\ 



the Cauchy- Schwartz inequality. The equal case ||M||i = ||M|| is 
derived from M'M = I or AP = M. 

Moreover, the main theorem is more efficient than previous 
results'^ Our result shows that the approximation error for syn- 
thesizing quantum qudit operation can be close to zero in terms of 
nonlocal quantum gate cost. However, it'-" '^ reads 
.D{I,U) 



\\U-Ua\\<- 



2" 



-2D(7,[/),5-'(e^-(l+^' 



+ C2D{I,U}n*d' (31) 



which is inexplicit because D{I, U) depends p and the sum of the last 
two terms may be lower bounded with nonzero constant by choosing 
5 = l/(n^D{I, U)). Therefore the approximation interval length d 
should be smaller such as 1/(«'D(I, LO) with k> 3. Moreover, their 
error bounds have not shown explicit relationship between the 
approximation error and total nonlocal gate cost. Our result presents 
that the number of approximate nonlocal quantum operations is the 
order of with the approximation error bound g. In comparison to 
the previous results'* the present geometric way shows the rela- 
tionship between the quantum approximation and the evolution of 
special geodesic equation. With the quantum circuit model the 
detailed circuit has to be found for generating U while one needs 
to find the shortest geodesic curve linking I and U for geometric 
method. So, these results may be considered as different quantum 
approximations with their own features. 

Methods 

The Lemmas 2-6 are used to prove Main Theorem. Let H{t) (at time [0, T]) be the 
time-dependent normahzed Hamiltonian generating the minimal geodesic of length 
d{I, U) — T. Let ^2(0 be the projected Hamiltonian with only the one- and two-body 
items and generating U2 with evolution time T. From the Lemma 3 it follows that 



\\U-U2\\<p-''^cl"T 



(32) 



Divide the time interval [0, T\ into N equal parts with length d = T/N. Let U2 be the 
unitary generated by H2{t) with evolution time S in the j-th time interval, and C/^ be 



(t)iitwith 



the unitary generated by the mean Hamiltonian H - 

evolution time d. F{H2{t)) = ^ ^ . hj < 1 because F{Hf = h- -Hp , h-j = 1 

with reasonable scaling. Hence ||H2(t)|| = ||y^^. ?ijffj|| < \hj\ <L(^y^^_^ h^'j 

<Sn where I — d^{(f — l){n^ — n)l2 is the number of one- and two-body items in 
H{t), Le. the number of terms in H2{t). Then from the Lemma 2 it follows that 



m-u' kc- 



(33) 



Moreover, from the Lemma 6 there exists a unitary L^^ synthesized using at most 
Nd^n^ one- and two-qudit gates, and satisfies 



\U>M-U>\ 



-'"'^-U>,\\<cn^T^N-' 



(34) 



with bounded on every time interval and c = cfn. 

In the follow we show that how to construct U2 and Uj^ in terms of and [/^ , 
respectively. In fact, U2 is generated using of H2 ( t) as 



dt 



Ui(t)^-iHi{t)Ui{t) 



(35) 



with Ul{0)^I,U^{{k-l)S)^U2~\^nd defined on [(j - 1)3, jd] for; - l/--,N. 
Letting be an identity operator on time interval [0, NS] except [{j — 1)6, jS]. 
(72-[/^C/^-^ •• • U2 is generated by the Hamiltonians • • ■ ,H^(f). L^^canbe 

generated similarly. Therefore 

\\U2-uj = \\u^uli-'- 



(36) 
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(37) 



From equations (33), (34) and (36) the approximation error 

||U-[/A||<||[/-t/2|| + ||[/2-f/A|| 
N 

<p-'^^d"T-\-Y,\\ui-U{\\ 

; = i 

<p-^nd"T^f2{\H-U^\ + \\ui,-u{\ 

^p-'^^d^T+TS-' {^e'''"-{l^d'n3)^ vein's- 
^p-^'^d^T + CaTd^n^d^Cin^b^ 



for a constant c - d"do + Codod'^n^ + c ^n^, where e'^^ "'^ - {l+d^nS) ^0{d'*n^S^) and 
Co is a constant. Choose a suitable penalty p so that D{I, U) = T satisfies Sdo/S ^ T 

T 

c/q from the Lemma 4 and]? ^'^ :^ £, and small ^ — ^ {N^ cc ) satisfying S e. Since 

C/^ may be synthesized with Nd^ri^ one- and two-body gates the total number is 
N^d^n^ = 0{e~^) one- and two-body gates. 
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